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Closed-fom  and  finite-element  solutions  are  presented  for  the 
themoelastic  behavior  of  laminated  composite  shells.  The  material 
of  each  layer  is  assumed  to  be  thermoelastioally  orthotropic  and  bi- 
modular,  i.e.,  having  different  properties  depending  upon  whether 
the  fiber-dimction  normal  strain  is  tensile  or  compressive.  The 
formulations  are  based  on  the  themoelastic  generalisation  of  Dong 
and  Tso ' 8  laminated  shell  theory ,  which  includes  thickness  shear 
deformations.  The  finite  element  used  here  has  five  degrees  of 
freedom  per  node  (three  displacements  and  two  bending  slopes). 
Numerical  results  are  presented  for  deflections  and  the  positions 
of  the  neutral  surfaces  associated  with  bending  along  both  coordi¬ 
nate  directions.  The  closed- form,  and  finite-element  results  are 
found  to  be  in  good  agreement. 
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INTRODUCTION 

As  the  field  of  composite-material  mechanics  becomes  more  highly  developed. 
Increasing  attention  Is  being  given  to  the  development  of  more  realistic 


models  of  actual  material  behavior  and  to  the  application  of  these  models  to 
thermostructural  analysis  of  composite-material  structural  elements  such  as 
plates  and  shells.  Certain  fiber- reinforced  composite  materials,  especially 
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those  with  very  soft  matrices,  exhibit  the  Interesting  phenomenon  of  having 
quite  different  elastic  properties  when  loaded  along  the  fiber  direction  In 
tension  as  contrasted  to  compression.  This  was  demonstrated  for  cord-rubber 
composites  by  Clark  [1]  and  Patel  et  al.  [2].  The  first  attempt  to  formulate 
a  theory  of  elastic  behavior  of  such  materials  was  due  to  Ambartsumyan  [3], 
and  a  comprehensive  theory  consistent  with  experimental  results  was  Introduced 
In  [4]  and  further  discussed  in  [5]. 

Most  of  the  thermoelastic  analyses  of  bimodulus-material  structural 
elements  [6-12]  have  been  limited  to  isotropic  bimodulus  materials.  However, 
in  a  recently  developed  theory  of  micromechanics  of  fiber-reinforced  materials 
with  soft  matrices  [13],  it  was  shown  that  the  thermal -expansion  coefficients, 
as  well  as  the  elastic  properties,  should  depend  upon  the  sign  of  the  fiber- 
direction  strain.  Unfortunately,  there  does  not  appear  to  be  any  appropriate 
experimental  data  available  to  date  to  confirm  this  conclusion  and  to  provide 
quantitative  values  for  the  thermal -expans ion  coefficients.  To  the  best  of 
the  knowledge  of  the  current  Investigators,  the  only  analysis  to  provide  for 


the  bimodular  effect  on  thermal  expansion  is  a  very  recent  thermoelastic 
analysis  of  thick  laminated  plates  [14]. 

There  have  been  a  few  thermoelastic  analyses  of  shells  laminated  of 
ordinary  composite  materials.  Stavsky  and  Smolash  [15]  considered  thin, 
laminated,  orthotropic  shells  using  Love's  first-approximation  shell  theory, 
and  Pao  [16]  treated  similar  shells  using  Fliigge's  higher-order  thln-shell 
theory.  Recently,  Padovan  and  Lestingl  [17]  analyzed  heated  anisotropic 
shells  Including  thickness-shear  deformation. 

A  number  of  analyses  of  various  kinds  of  bimodulus  shells  have  appeared 
In  the  literature;  ten  of  them  were  reviewed  In  [18].  However,  In  every 
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Instance  they  treated  only  thin  shells  subjected  to  mechanical  loading  only. 

The  present  analysis  Is  believed  to  be  the  first  analysis  of  blmodulus 
shells  to  Include  either  thermal  loading  or  thickness-shear  deformation. 

The  theory  used  Is  a  generalized  first-approximation  thermoelastlc  shell 
theory  which  can  be  reduced  by  means  of  tracer  coefficients  to  various 
simpler  theories. 

GOVERNING  EQUATIONS 

Let  x  and  y  denote  the  axial  and  circumferential  position  coordinates 
measured  on  the  shell  middle  surface,  and  z  the  outward  normal  position 
coordinate  (see  Figure  1). 

The  displacement  field  at  an  arbitrary  location  (x,y,z)  is  given  by 
U(x,y,z)  *  u(x,y)  +  zex(x,y) 

V(x,y,z)  *  v(x,y)  +  z By(x,y)  (1) 

W(x,y,z)  *  w(x,y) 

Here,  u,v,w  are  the  middle- surface  displacements,  and  a  and  a  are  the 

A  Jr 

bending  slopes. 

The  strain-displacement  relations  for  small  deflections  can  be  written 
as 

ei  -  +  zKi  (1*1 ,2, 4, 5, 6)  (2) 

Here,  are  the  engineering-strain  components  at  an  arbitrary  location 
(x,y,z),  ej  are  the  engineering -strain  components  on  the  middle  surface 
(x,y,0),  and  kj  are  the  curvature  changes.  The  notation  of  classical  com¬ 
posite-material  mechanics  Is  used,  with  1  and  2  denoting  normal  action  In 


directions  x  and  y,  respectively,  and  6  denoting  shear  action  with  respect 
to  x.y  axes.  Now 


e$  *  u’x  *  e2  *  V*y  +  WR)  •  *  U»y  +  v*x 

e4  3  6y+w*y"  (c1/R>v  *  e°  55  6x  +  w’x 

*1  3  6x,x  ’  <2  *  ey,y  *  *6  *  ex,y  +  ey,x 

k4  *  k5  =  0  +  (C2/2R)(v.x-u,y) 


Here,  R  Is  the  radius  of  the  middle-surface,  the  are  shell-theory  tracers 
to  be  discussed  later,  and  (  ),x  s  3(  )/3x. 

Considering  the  shell  to  consist  of  either  a  single  orthotropic  layer 
or  to  be  a  cross-ply  laminate  (one  having  all  layers  oriented  at  either  0° 
or  90°  with  respect  to  the  cylinder  axis),  the  thermoelastic  version  of 
generalized  Hooke's  law  may  be  written  as  follows  for  each  layer: 
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Here,  are  the  stress  components,  are  the  plane-stress-reduced 

stiffnesses,  oj^  are  the  thermal  expansion  coefficients,  and  T  Is  the 
temperature  measured  from  the  strain-free  temperature.  Subscript  k*l  for 
fiber-direction  tension  and  2  for  compression,  and  subscript  l  denotes  the 
layer  number. 

The  stress  resultants  and  stress  couples  are  defined  In  the  customary 
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May  for  a  first-approximation  shell  theory  as 

h/2 

(N1,M1)  -  (l.z)^  dz 

-h/2 


(5) 


where  h  =  total  laminate  thickness.  Similarly,  thermoelastic  stress  resul¬ 


tants  and  stress  couples  are  defined  as 

h/2 

-  f  (l,z)E1fcedz 
-h/2 


where 


Em  *  Qijwajw  (no  sum  on  u) 


(6) 


(7) 


Since  thickness-shear  deformation  Is  included,  the  shear  stress  re¬ 
sultants  are  introduced 


(Qz»Ql)  *  I  (crv»®5)dZ 

-h/2 


(8) 


Substituting  Equations  (4)  and  (6)  into  Equations  (5)  and  (8),  one 
obtains  the  following  shell  constitutive  relations 
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Here,  ■  Inplane  stiffness,  B^  *  Inplane-bending  coupling  stiffness, 
Djj  *  bending  stiffness,  S.^  *  thickness  shear  stiffness,  defined  by 

h/2 

{AirB1j»D1j)  *  J  (l.z.z2)Qij  dz  (1 ,4*1 ,2.6) 


-h/2 
h/2 


(10) 


S1j  "  K2J 
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(1  J-4,5) 


(Derivations  of  A^,  B^,  D^,  nJ,  and  are  carried  out  In  detail  for  a 
laminated  bimodulus  shell  In  Appendix  A.) 

The  shell  equilibrium  equations.  In  the  absence  of  body  forces  and 
body  moments,  can  be  written  as 
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where  P  is  the  normal  pressure. 

If  the  shell-theory  tracer  Ci  is  set  equal  to  unity  and  C2  set  equal  to 
zero,  the  theory  presented  here  can  be  considered  to  be  the  thermoelastic 
version  of  the  shear-deformation  shell  theory  of  laminated,  orthotropic, 
circular  cylindrical  shells  presented  by  Dong  and  Tso  [19].  This  theory  Is 
the  shear-deformable,  laminated,  orthotropic  version  of  the  well-known  Love 
first-approximation  shell  theory  [20],  as  modified  by  Reissner  [21];  see 


If  the  shell-theory  tracers  are  set  equal  to  other  values  as  specified 
In  Table  1,  the  theory  represents  the  shear-deformable,  laminated  orthotropic 
version  of  the  Sanders  "best"  first-approximation  theory  [23],  Loo's  approxi¬ 
mate  theory  [24],  Morley's  shallow-shell  theory  [25],  and  Donnell's  very- 
shallow-shell  theory  [26].  It  Is  Interesting  to  note  that  when  generalized 
to  Include  shear  deformation,  one  cannot  distinguish  between  Love's  first- 
approximation  theory  and  Loo's  theory  and  also  between  Morley's  and  Donnell's 
shallow-shell  theories. 

Substituting  Equations  (3)  and  (9)  into  Equations  (11),  we  obtain  the 
following  operator  equation 

[L]{5}  -  {f>  (12) 

where 

(d>  a  {U,V,W,&x,ey}T 

Now  [L]  Is  the  symmetric  matrix  of  the  following  differential  operators: 

L11  *  Alldx  +  (A66  ■  ^2^66  +  |  ^2®66)dy 
Li2  *  (A12  +  A66  C2D66)dxdy 
Li3  ■  (A12/R)dx  ;  L14  *  Blldx  +  (®66  ^2B66)dy 

Lis  *  (B12  +  B66  -  y  C2D66)dxdy 

L22  ■  (A$$  +  C2B66  S2D6g)dx  +  ^22dy  "  ^XB44 
L23  *  (R  XA22  +  CiS44)dy  ;  L24  *  (B^Bee*^  C2066)dxdy 

L2s  •  (Bgs+'j'  C2066)dx  +  B22dy  +  CjSm, 

L33  *  -  Sssd*  -  S44dy  +  (A22/R  ) 

L31,  *  (R  lBi2"Sss)dx  ;  L35  *  (R  1B22*S44)dy 


(13) 
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L44  ■  Dud*  +  066d*  -  S55  ;  L45  a  (D12  +  066)dxdy 


1*55  *  ^66<lx  +  ^22^y  "  S44 


£,  -  C^/R 


dx  *  3(  )/3x,  etc. 


Also  the  components  of  the  generalized  thermal-force  vector  {f}  are: 


fi“Ni,x  1 


f 2  -  N^y  i  f3  •  P  -  (N^/R) 

f  *  mT  ;  f  *  MT 
4  i,x  5  2 ,y 


(14) 


In  view  of  the  assumed  linearity  of  the  displacements  with  z,  It  is 
consistent  to  assume  that  the  temperature  distribution  is  also  linear  with  z: 


T(x,y,z)  •  T  (x,y)  +  zT^x.y) 


(15) 


CRITERIA  FOR  HOMOGENEITY  ALONG  MIDDLE  SURFACE 

In  deriving  Equations  (12),  we  tacitly  assumed  that  the  laminate  stiff¬ 
nesses  (Aij»Bij»Dij»Sfj)  are  *11  independent  of  coordinates  (x,y)  on  the 

middle  surface.  However,  in  view  of  the  bimodulus  nature  of  the  materials 

• 

comprising  the  laminate,  these  stiffnesses  depend  upon  the  fiber-direction 

neutral -surface  positions  associated  with  the  respective  layers  (i.e.,  znx 

for  a  single  layer  with  axially  oriented  fibers,  and  z„w  and  z_„  for  a  cross- 

nx  ny 

ply  laminate). 

Thus,  for  layers  having  the  fibers  oriented  axially,  the  associated 
fiber-direction  neutral -surface  position  is  determined  by 


£i  *  ei  + 


or 


nx 


-  Cj/k!  a  -  u*x/eXfX  *  constant 


(16) 


Similarly,  for  layers  having  the  fibers  oriented  circumferentially 

62  ■  4  +  znv*2  “  0 


or 


zny  *  ‘  e2/K2  ’  (v,y  +  R"lw)/By.y  *  constant 


(17) 


CLOSED-FORM  SOLUTION 

A  solution  is  sought  which  satisfies  the  governing  operator  equation. 
Equation  (12),  the  subsidiary  relations.  Equations  (16)  and  (17),  and  the 
appropriate  boundary  conditions. 

A  closed-form  solution  has  been  found  for  the  following  conditions: 
Loadi ng  (sinusoidally  distributed): 

P  *  PQ  sin  ax  sin  By  ,  TQ  *  TQ  sin  ax  sin  By 

Tx  =  ij  sin  ax  sin  By  »  a  *  irx/a  ,  B  *  iry/b 


Here  a  and  b  are  the  dimensions  of  the  shell  In  the  respective  axial  and 
circumferential  directions  (see  Fig.  1). 

Boundary  Conditions  (freely  supported): 

Ni(0,y)  -  Nj(a,y)  *  M^O.y)  «  Mx(a,y)  *  0 

w  (0,y)  =  w(a,y)  «  v(0,y)  =  v(a,y)  =  0 

N2(x,0)  =  N2(x»b)  =  M2(x,0)  =  M2(x,b)  «  0  (19) 

w  (x,0)  *  w(x,b)  *  u(x,0)  »  u(x,b)  *  0 

By(0,y)  *  sy(a,y)  *  Bx(x,0)  =  Bx(x,b)  =  0 

Under  these  conditions,  the  solution  to  Equation  (12)  Is  of  the  form 

u(x,y)  *  0  cos  ax  sin  By 

v(x,y)  *  V  sin  ax  cos  By 

w(x,y)  *  W  sin  ax  sin  sy 

Bx(x,y)  ■  X  cos  ax  sin  sy 

By(x,y)  *  Y  sin  ax  cos  sy 


(20) 
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Substitution  of  Equations  (20)  Into  Equation  (12)  leads  to  the  follow¬ 
ing  nonhomogeneous  algebraic  system: 


[C]U}  *  CF} 


where 


{A}  »  {U,V,W,X,Y> 

{F}  *  {F1,F2,F3,F4,F5}T 

The  quantities  Fr  and  the  coefficients  CfS  of  the  matrix  [C]  are  not 
presented  here,  for  brevity. 

For  a  given  set  a,  6,  PQ,  R,  F..  and  either  single-layer  or  cross-ply 
construction,  one  needs  to  solve  the  5x5  matrix  Equation  (21)  for  the  vector 
{a}  of  amplitudes  of  the  generalized  displacements,  subject  to  subsidiary 
conditions  (16)  and  (17).  For  bimodulus-material  shells,  the  laminate  stiff¬ 
nesses  (A.j,  Bij*Dij)  are,  in  general,  not  constant,  but  depend  upon  x  and  y 
through  the  fiber-direction  neutral-surface  positions  (znx  and  zny).  How¬ 
ever,  for  the  present  combination  of  loading  and  boundary  conditions,  znx 
and  zny  are  both  constants,  i.e.,  independent  of  x  and  y.  Although  it  is 
conceptually  possible  to  substitute  the  solution  functions  into  Equations 

(16)  and  (17)  to  obtain  cubic  equations  involving  z  and  z  ,  it  is  computa- 

••a  ny 

tionally  much  more  efficient  to  satisfy  Equations  (16)  and  (17)  by  iterating 
on  zrtx  snd  V 


FINITE-ELEMENT  FORMULATION 

Since  an  exact  closed-form  solution  to  Equation  (12)  can  be  obtained 
only  under  special  conditions  of  geometry,  edge  conditions,  loadings,  and 
lamination,  it  is  desirable  to  have  available  a  more  general  method.  Here, 
we  develop  a  simple,  mixed-type,  finite-element  formulation  which  has  no 
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such  limitations,  except  for  those  Implied  In  the  formulation  of  shear- 
flexible  laminated  shell  theory. 

Let  the  region  R  be  subdivided  Into  a  finite  number  N  of  subregions: 
finite  elements,  Rg  (e*l,2,...,N).  Over  each  element,  the  generalized 
displacements  (u,v,w,s  ,e  )  are  Interpolated  according  to 

X  jr 

r  1  r  1  S  2 
U  *  £  U^  ,  V  *  £  Vi^i  ,  W  ■  £  W^. 

(23) 
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Here,$“  (o=1,2,3)  Is  the  Interpolation  function  corresponding  to  the  1-th 
node  In  the  element.  It  Is  noted  that  the  Inplane  displacements,  the 
normal  deflection,  and  the  bending  slopes  may  be  approximated  by  different 
sets  of  Interpolation  functions.  Although  this  generality  Is  considered 
in  the  formulation  presented  here,  when  the  element  is  actually  progranmed, 
we  set  ^  1»^2=<|> 3  (r*s»p)  for  simplicity.  Noting  that  r,  s,  and  p  denote 
the  number  of  degrees  of  freedom  (OOF)  for  each  variable,  the  total  number 
of  DOF  per  element  Is  2r  +  s  +  2p. 

Substituting  interpolations  of  the  form  (23)  for  u,  v,  w,  sx,  and 
Into  the  Galerkln  Integrals  associated  with  the  governing  operator  equation 
(12),  we  obtain 

([L]{5>  -  {f»{*)dxdy  »  0  (24) 

Now  using  Integration  by  parts  once  in  order  to  distribute  the  differentiation 
equally  among  the  terms  In  each  expression,  we  obtain  the  element  equation 

CK]e{We  -  {P>e  (25) 


The  elements  K®?  (a,6a1 ,2,. . . ,5)  of  the  element  stiffness  matrix  and 
of  the  generalized  force  vector  are  listed  In  Appendix  B. 

In  the  present  Investigation,  cylindrlcally  curved  rectangular  elements 
of  the  serendipity  family  are  used,  with  the  same  interpolation  for  all  of 
the  variables.  The  resulting  stiffness  matrices  are  20  by  20  for  the  four- 
node  element  and  40  by  40  for  the  eight-node  element.  Reduced  integration 
[27,28]  must  be  used  to  evaluate  the  matrix  coefficients  in  Appendix  B. 

For  example,  for  the  four-node  rectangular  element,  the  1x1  Gauss  rule  must 
be  used  rather  than  the  standard  2x2  Gauss  rule. 

NUMERICAL  RESULTS 

It  would  be  desirable  to  compare  the  results  obtained  by  the  present 
analyses  with  those  given  in  the  literature  for  special  cases.  Unfortunately, 
however,  there  is  a  dearth  of  solutions  of  cylindrical  panels  subjected  to 
sinusoidally  distributed  mechanical  and  thermal  loadings.  However,  in  a 
recent  closed-form  and  finite-element  study  [14]  of  thermally  loaded  plates, 
good  agreement  was  obtained  with  results  presented  by  Boley  and  Weiner  [29] 
for  isotropic,  thin  plates. 

As  practical  examples  of  orthotropic  bimodulus  materials,  the  same 
two  unidirectional  cord-rubber  materials  as  considered  in  [14,18]  are  con¬ 
sidered,  namely,  aramid-rubber  and  polyester-rubber.  The  inplane  elastic 
properties  were  obtained  from  experimental  results  of  [2]  using  the  data- 
reduction  procedure  presented  in  [4].  Since  the  thickness-shear  moduli  were 
not  measured  in  [2],  they  were  estimated  as  described  in  detail  in  [30]. 

The  elastic  properties  are  listed  in  Table  2. 

Unfortunately,  the  present  investigators  are  unaware  of  any  experimentally 


determined  values  for  the  thermal -expansion  coefficients  of  cord-rubber 
materials.  However,  the  mlcromechanlcs  analysis  of  blmodular  action  pre¬ 
sented  In  [13]  suggested  that  the  thermal -expans Ion  coefficients  (ax  and  a2) 
of  these  materials  should  also  depend  upon  the  sign  of  the  fiber-direction 
strain.  Thus,  In  the  numerical  calculations  presented  here,  the  following 
dimensionless  relationships  are  used: 

ax/ai  *  0.5  ;  a2/a§  «  1.0  ;  af/a*  »  0.1 

Here,  superscripts  c  and  t  refer  to  compressive  and  tensile  fiber-direction 
strains. 

Table  3  shows  the  effect  of  the  radius-to-thickness  ratio  (R/h)  on 
the  locations  of  neutral  surfaces  and  dimensionless  deflections  for  single¬ 
layer  and  two-layer  cross-ply  aramld-rubber  cylindrical  panels  under  sinu¬ 
soidal  mechanical  loading  by  Sanders  theory.  As  the  radlus-to-thlckness 
ratio  Is  Increased  to  Infinity,  the  panel  can  be  considered  as  a  plate. 

Table  3  also  shows  the  convergence  of  the  dimensionless  deflections. 

Numerical  results  of  the  Influence  of  the  aspect  ratio  on  the  dimen¬ 
sionless  deflections  and  neutral -surface  locations  for  single-layer  and 
two-layer  cross-ply,  freely  supported  cylindrical  shells  constructed  of 
bimodulus  materials  and  subjected  to  sinusoidal  thermal  loading  (R/h  •  10) 
by  Sanders  theory  are  shown  In  Tables  4  and  5,  respectively.  Again,  there 
Is  a  close  agreement  between  the  finite-element  and  closed-form  results. 

Figure  2  shows  the  effect  of  radlus-to-thlckness  ratio  arid  aspect 
ratio  (a/b)  on  the  dimensionless  deflections  and  neutral -surface  locations 
for  one-layer  and  two-layer  cross-ply,  freely  supported  aramld-rubber 
cylindrical  shells  under  sinusoidal  thermal  loading  by  Sanders  theory. 
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Figures  3  and  4  show  the  Influence  of  aspect  ratio  and  radlus-to- 
thlckness  ratio,  respectively,  on  the  locations  of  neutral  surfaces  for 
single-layer,  freely  supported,  aramld-rubber  cylindrical  shells  under 
sinusoidal  thermal  loading  by  Sanders  theory.  Similar  results  are  shown 
In  Figures  5  and  6  for  two-layer  cross-ply,  freely-supported  polyester- 
rubber  cylindrical  shells  under  sinusoidal  thermal  loading.  As  can  be 
seen,  there  Is  only  a  slight  change  of  neutral -surface  locations  for 
radlus-to-thlckness  ratios  greater  than  60. 


CONCLUDING  REMARKS 


' 


Jv  - 


A  finite-element  formulation  of  equations  governing  layered  anisotropic 
composite  shells  subjected  to  mechanical  as  well  as  thermal  loading  is 
presented.  The  element  includes  the  effect  of  shear  deformation  and  in¬ 
volves  five  degree  of  freedom  (three  deflections  and  two  rotation  functions) 
per  node.  Numerical  convergence  of  linear  and  quadratic  elements  Is  shown, 
and  results  are  presented  for  single-layer  and  two-layer  cylindrical ly  curved 
cross-ply  panels  subjected  to  sinusoidal  and  uniform  loadings:  thermal, 
mechanical,  and  combined  loadings  are  considered. 

To  check  the  finite-element  results,  a  closed-form  solution  is  developed 
herein  for  cross-ply  cyllndrlcally  curved  panels  subjected  to  sinusoidal 
mechanical  and/or  thermal  loadings.  The  exact  solution  can  be  obtained  only 
under  special  conditions  of  geometry,  edge  conditions,  and  loadings.  How¬ 
ever,  the  finite-element  formulation  presented  here  does  not  have  any  limi¬ 
tations  except  for  those  Implied  In  the  formulation  of  the  governing  equations. 
The  finite-element  solutions  are  found  to  be  In  close  agreement  with  the 
closed-form  solutions  for  2  by  2  mesh  of  quadratic  elements  in  the  quarter 
shell . 
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Thus,  the  finite  element  developed  here  Is  computationally  simply 
compared  to  other  cylindrical  shell  elements  used  previously  In  the  thermal 
stress  analysis  of  cylindrical  shells. 

Extension  of  the  present  element  to  non-linear  thermal  stress  analysis 
and  to  thermal  buckling  analysis  Is  recommended.  In  those  cases,  the 
present  element  should  result  In  substantial  savings. 
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APPENDIX  A 

DERIVATION  OF  EXPRESSIONS  FOR  THERMAL  FORCES  AND  MOMENTS 


Cast  I 

For  Case  I.  >  0  and  zfty  <  0  with  zM  governing  loyer  1  (O*)  and  z^ 
layer  2  (90*). 

T  f2W 

"x  *  1^2^1122  ®122  +  <J1222  °222)T  d2  +  Jz  *QU12  a112  +  <l1212  °212)T  dZ 

rznx  .  fh/Z 

♦  JQ  {Q1121  “121+Q1221  “221)T  d2  +  J  (Qmi  a111+Q1211  “211)T  dz  (A"1) 

u  znx 

Let 

(Q1122  °122  +  Q1222  °222}  *  B122  *  (()in2  °112  +  Q1212  a212)  "  B112 

(A-2) 

*Q1121  a121+Q1221  “221  *  "  B121  *  *Q1111  “111+Q1211  “211  *  *  Blll  •  etc* 

Then, 

Nx  *  [b122  V2ny  +  h/2)  +  ®112  V°‘2ny*  +  S121  V2nx~0) 

+  Bin  Vh/2-z^)  ♦  B122(T1/2H)(z^.h2/4) 

+  WV2^0'2^  +  812l(T1/2h)(z5x-0) 

♦  BmtT^hHh^A-z^JlsIn  ax  sin  By 

Nx  “  ^(®122 +  8m)(T0h/2^  +  (B121"Blll)  Vnx  *  ^122**112*  To2ny 
+  {®m-®i22>^Tih/8>  ♦  <Bi2i-Bni)<Ti2ix/2h) 

♦  (*122’8n2^Tl2ny/2h^Sln  **  ,1n  By  (A-3) 


Similarly, 


20 


i 

i 


■  E(*222  *  8211^Toh^  +  ^®221  “  8211^  *oznx  +  ^®222“8212^  Tozny 

+  ^®211  "  ®222^1^8^  +  ^®221 "  ®211^1znx^2^  +  ^®222”®212^ 
(T1zjy/2h)]s1n  ax  sin  By 


(A-4) 


T  fh/2  f0  nx  h/2 

Mi  *  J_h/28122  Tz  d2  +  J2  6112Tzd2  +  J0  8121Tzdz+  2  8niTzdz 

ny  nx 

[(Bin  -  B122)(T0hZ/8)  ♦  ($121  -  ®i , -| ) (T02S*/2)  +  <8122-  8112)(ToZV2) 

+  ^8122  +  8m^T1h2/24)  +  ^®121  *  ®111^1zru/3^ 

♦  (*122~8n2,(Tlzny/3h^sin  “*  s1n  8y 
Similarly, 

M J  *  C(e21re222)(Toh2/8J  +  <822i-62n)(To4/2>  + 

♦  (8222  +  ®211  ^^l^1*^24^  +  ^ 8221  ”  ®2]]^^lZnx^3^ 

♦  (B^- 8212)(T1zJy/3h)]s1n  ax  cos  By 

Using  the  above  equations  In  conjunction  with  equations  (12)  and  (18),  we 


(A— 5) 


(A-6) 


obtain  the  following: 


Rx,x  *  a^®122  +  ®111^Toh/2*  +  *®121 "  ftm^Toznx  +  *8122' ®112*Vr 


+  (®m  *  8122^1h^  ♦  (8721  "  ®m^^1znx^2h^ 
♦  (B122-  8n2^Tlzny/2h^ 


(A-7) 


% 


(A-8) 


+  {®211  *  B222^^1^8^  *  (^221  *  B211^1zn?/2>>^ 

♦  (S222  ‘  B212^Tlzn/2h^ 

**x,x  *  «C<8ni  *  ®122^Toh2/8)  +  (B121’Blll^Toznx/2*  +  *B122“6112* 

^Vo/2*  +  *b122  +  Blll){Tlh2/24)  +  (Bi2l*Bni^Tlznx/3h^ 

+  ^b122  “  B112^1zny^3^  (A"9) 

Ry,y  *  B^B211  ■  B222^Toh2/8)  *  *b221  '  B211^Toznx/2)  *  ta222  “  B212^V5/2) 

*  (®222+  b211^^1^  /^4)  +  ( B221  *  B211^^1znx^3^  +  ^B222”  B212^^1zny^3^^ 

( A— 1 0) 

In  a  similar  way,  one  can  obtain  the  expressions  for  the  above-mentioned 

quantities  for  the  remaining  seven  cases  as  follows: 

Case  II  (z^O,  zny>0) 

*x,x  *  oC(b122  +  Blll^Toh/2)  +  ^®121  *Blll)(Vnx) 

+  (ain  -  8122)(T1h/8)  +  (0121  -  Bm^Tlznx/2h)l 

fiy,y  *  b^b222  +  B211)(Toh/2)  +  ^ *221  "  B211^ToZnx^ 

+  <b211  *  B222^*lh/8)  +  *B221  '  B211^T1znx/2h^  (A-11) 

*x,x  *  *^Blll’B122^Toh2/8^  +  *®121  “Blll^Vnx/2^ 

+  t®122  +  Blll ) 1  ®**/24)  ♦  (0^21  “  ®lll )(^lznx/3h)^ 

fiy,y  ’  b^b211  *  ®222^Toh2/8)  *  *b221  "  B211^Toznx/2) 

+  (®222"  b211  )(^i^2/2^)  +  ( b221  *  B211^1znx^3^ 
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Case  III  (z^O,  z^O) 

C  "  °^B122  +  Blll)(Toh/2)  +  ^121  ’e111^Vny^ 

♦  (Bj22“  B112^oznx^  +  ^111  “  s122^^1h/8^ 

+  (ftl 21  *®m^Tizny/2h)  *  *Bl22’B112^Tlznx/2h^ 

^y»y  *  B^B222  +  +  ^221  ‘B211^ozny^ 

+  (B222*B212^Vnx)  +  *B211  "  B222)(Tlh/8) 

+  ^b221  “  B211^1zny^2h^  +  ^B222  “  8212^1ziu/2h^ 

RI,x  *  «t(Blir8122)(t0h2/8)  ♦  (812l-81n)(Tozny/2) 

♦  (8122*B112^Toznx/2)  +  *B122  *  B111^1hJ/24^ 

+  (B12l'Blll)^12J/3h)  +  (B122“B112^Vnx/3h^ 

^y.y  *  ^b211  ”  b222^V*  ^  +  ^  b221  ■  82n)(ToZn/2) 

♦  U222-  B212^Toznx/2*  +  ^B222  +  B211^T1h  /24* 

+  ^B221  "  B211^^1zn/3h^  +  *B222"  B212^1Znx/3h^ 

^SSJl  <znx<0’  zny<0) 

C  *  ®^b122  +  B111)(Toh/2)  +  ^b122"  B121^ozny^ 

+  <8111  “  B122)(Tlh/8)  +  ^8122  *  B121^^1zny^2^^ 


RJ,y  “  8C^®222  *  *211  +  ^B222”  B221^ozny^ 

♦  {®211 *  B222^1^8^  +  ^b222”  B221^1zny^2^ 


CA- 13) 


23 


•  y  ■ 


•f-, , 

'  V  ' 


8?. 


fe 


lifts 


»•"* 


i 


w 


'%■ 


1 


ftx,x  *  *^8m  “  8122^Toh*/8)  +  *8122‘  8121^V5/2) 

+  ^8122  *  8in)<Tlh  +  ^8122”  8121^1zpy^2^ 

ftJ,y  *  8^8211  “  8222^o^8^  +  ^8222~  8221^ozny^ 

*  (#222  +  8211 )  (^1^/24)  +  (®222*  ®221^^1*n/^^ 
For  neutral  surface  going  out  of  plane. 


Case 


V  (W0.5.  2^-0. 5) 


x,x 


N 


y.y 

lT 

x,x 


y.y 


*  °t(Si21  *  8112^^/2^  +  ^8121  “ 

*  SC(  &221  +  8212^^/2^  +  ^8221  * 

*  <*[Ui2i  “  *112^^0^®^  +  ^8121  + 
■  0[(0221  *  ®212^^/8^  +  ^8221  + 


8112^ty8^ 

8212^1^8^ 

8112)(Tl/24)] 

8212)(Tl/24)] 


Case 


VI  (2nx<“°*8*  V°*5) 


»i.x 

"J.y 

aT 

"x.x 

rt 

y.y 


*  <*[(0121  +  8112^^/2^  +  ^8121  " 

*  0[(022l  +  8212^o^  *  ^8221  " 

*  <»[(0i2i  *  8112^^/8^  +  ^8121  + 

*  0[(®221 "  8212^^/8^  +  ^8221  + 


8112)(V8>] 

8212^1^8^ 

8ii2)(Ti/24)] 

8212)(ti/24)3 


Case 


VII  (z^O.5,  zny>0.5) 


x,x 


y.y 


x,x 


Rt 

y.y 


■  *[(®m  +  8h2^V2)  +  *8iu ' 

•  0[(02n  +  8212)(To/2)  +  {82U“ 

■  ®[(®in  *  Bii2)(To/8)  +  (8m  + 

■  0[(02n  “  0212^V8^  +  ^®211  * 


8n2)(Ti/8)3 

82i2^(^i^8^3 

8112,(V24)] 

8212)(T1/24)] 


ilL 


k; 
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Case  VIII  (znx<-0.5,  ^<*0*5) 

Rx,x  *  °^8121  +  8122^V2*  +  *8121  “  B122^V8^ 

Ny.y  “  bC(b221  +  8222)(To/2)  *  (b221  “  B222)(V8)] 

»I,x  -  «C(e121-8122)<y8>  +  (8121  *  6122)(T ^24)3 

»J#y  *  8C(8221  -  8222)(T0/8)  ♦  U221  ♦  6222)^/24)] 

For  a  single  layer,  change  8^2  to  8^*  8122  to  B12]»  6212  to  8211  and 
8 222  t0  6221 * 


(A-17) 
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APPENDIX  B 

LISTING  OF  COEFFICIENTS  OF  STIFFNESS  MATRIX 
AND  FORCE  VECTOR  FOR  FINITE-ELEMENT  FORMULATION 


The  elements  of  the  stiffness-matrix  appearing  In  Equation  (25)  are: 
Kij  "  AllG1j  +  *A66  "  C2B66+T  e2D66JG1j 

Kij  3  <VRKj  ;  Kij  3  BnHi,j  +  (B66  *  7  C2D66)H1j 

Kij  B12Hij  +  ^B66  ‘  7  C2D66)Hji 

Kij  3  A22Gij  +  ^A66  +  C2B66  +  \  C2D66^Gi j  +  G1S44Gij 

K1j  =  (A22/R)Mji  ;  Kij  3  (B66  +  7  C2D66)Hij  +  B12Hji 

K25  *  <B66  +  1  C2D66}Hi  j  +  *  GlS44Hij 

K33  *  S55Sij  +  S44Sij  +  (a22/r2)S1 j 

K34  3  S55R1j  +  ^B12/R)Rji  5  *35  3  +  (B22/R^Rji 

K44  “  °llTiJ  +  D66Tij  +  D66Tij  +  S55Tij  ;  *45  3  D12T1j  +  D66Tji 

K55  3  D66Ti j  +  D22Tij  +  S44Tij 

The  generalized-force  elements  appearing  In  Equation  (25)  are: 

f!  ■  L  <Nl,x  *  "ly  ■  1  c2<y>*!  dx 

*e 

i  ’  L  (N2.y  *  "six  *  7  c2"6.x>*]  dx  « 
re 

F1  *  |R  (P  "  N2/R)*  1  dx  dy 

•1  -  L  <"l  ,x  *  M6,y>*?  dx  •  F1  ■  L  <M2,y  *  MLx>*1  dx  dy 


(B-l ) 


(B-2) 


where 


II 

c-n 

UT»r- 

♦1  a\  n  d*  dy 

«  1 »£  j  .n 

Ke 

( i »jal >2,... >n) 

3: 

c~>.  Z3 

M 

a*  dx  dy 
-n  ’1.5  J.n 

Re 

( i*l ,2, . . . »r  ;  j*l .2, 

a 

’  *i.5*lnd’“* 

(i»l ,2,. . . ,r  ;  jal .2, 

s5n  , 

ij 

.  d*  dy 

Ke 

(i,j*l  ,2, . . . .s) 

L  ♦t.s’j.n  d*  dy 

*e 

(1*1,2*... >s  ;  j*l »2, 

13  ■ 

a?  .$]  dx  dy 

Ke 

(1J*l,2,...»s) 

(5,n*0,x,y) 

and  Gx*  *  G*.,  etc.  In  the  special  case  in  which  ♦{  * 
ij  ij 

matrices  in  Equation  (8-3)  coincide. 


Figure  1.  Shell  geometry 
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Figure  2.  Transverse  deflection  vs.  aspect  ratio  and  radius-to-thickness 
ratio  for  single-layer  and  two-layer  cross-ply  cylindrical 
panels  under  sinusoidal  thermal  loading  by  Sanders  theory 
(Material:  aramid-rubber) . 
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Figure  3.  Neutral -surface  location  vs.  aspect  ratio  for  single-layer 
cylindrical  shells  under  sinusoidal  thermal  loading  by 
Sanders  theory  (Material:  aramid-rubber,  R/h»10) 
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Figure 


Neutral -surface  location  vs.  radlus-to-thlckness  ratio 
for  single-layer  cylindrical  shell  under  sinusoidal 
thermal  loading  by  Sanders  theory  (Material : aramid- 
rubber,  a/b«l .0) . 
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Figure  5.  Neutral -surface  location  vs.  aspect  ratio  for  two-layer 
cross-ply(0°/90°)cy1indr1cal  shells  under  sinusoidal 
thermal  loading  by  Sanders  theory  (Material:  polyester- 
rubber,  R/h»10j. 


Figure  6.  Neutral -surface  location  vs.  radius -to -thickness  ratio 

for  two-layer(0°/90B)  cylindrical  shells  under  sinusoidal 
thermal  loading  by  Sanders  theory  (Material :po1yester- 
rubber,  a/b-1 .0)  . 
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Table  1.  List  of  Shell -Theory  Tracers  and  Their  Values 


Theory  (Thln-Shell  Theory  Generalized 
to  Shear-Flexible  Theory) 

C1 

C2 

Sanders 1 

1 

1 

Love's  first  approximation 
and  Loo's 

1 

0 

Morley's  and  Donnell's 

0 

0 

Table  2.  Elastic  Properties  for  Two  Tire-Cord/Rubber.  Unidirectional,  Blnddulus 
Composite  Materials* 


Aramld-Rubber 

Polyester-Rubber 

Property  and  Units 

k-1 

k-2 

k-1 

k-2 

Longitudinal  Young's  modulus,  GPa 

3. SB 

0.0120 

0.617 

0.0369 

Transverse  Young's  modulus,  GPa 

0.00909 

0.0120 

0.00800 

0.0106 

Major  Poisson's  ratio,  dimensionless* 

0.416 

0.205 

0.475 

0.185 

Longitudinal -transverse  shear  modulus,  GPa0 

0.00370 

0.00370 

0.00262 

0.00267 

Transverse-thickness  shear  modulus,  GPa 

0.00290 

0.00499 

0.00233 

0.00475 

*F1ber-d1  recti  on  tension  Is  denoted  by  k-1,  and  fiber-direction  compression 
6by  k-2. 

It  Is  assuaed  that  the  minor  Poisson's  ratio  is  given  by  the  reciprocal 
relation. 

It  Is  assumed  that  the  longitudinal-thickness  shear  modulus  Is  equal  to  this 


Table  3.  Effect  of  the  radlus-to-thlckness  ratio  (R/h)  on  the 
locations  of  neutral  surfaces  and  deflections  for 
single-  and  two-layer, cross-ply ,a rami d -rubber  cylindrical 
panels  under  sinusoidal  loading  by  the  Sanders  theory 
(Tj  ■  T  *  0,  a/b  ■  1,  b/h  ■  10,  material  I). 


R/h 

Layers 

Source 

w 

Zx 

z. 

R/h-*-  • 

1 

CF 

0.02094 

0.44205 

-0.16185 

(plate) 

FEM 

0.02093 

0.44205 

-0.1616 

2 

CF 

0.01982 

0.4384 

-0.03418 

FEM 

0.01981 

0.4384 

-0.03416 

1 

CF 

0.020246 

0.4408 

-0.1840 

100 

FEM 

0.020234 

0.4408 

-0.1838 

2 

■H 

0.01892 

0.43666 

-0.036686 

MESH 

0.01891 

0.43661 

-0.03666 

1 

CF 

0.01943 

50 

FEM 

0.01943 

EE9 

USSR 

2 

CF 

0.01793 

-0.03909 

FEM 

0.01793 

-0.03905 

1 

HB 

0.01900 

40 

0.01899 

2 

CF 

0.01742 

0.4341 

-0.04027 

FEM 

0.01741 

0.4341 

-0.04022 

20 

1 

CF 

0.01663 

I 

FEM 

0.01663 

2 

mmm 

0.01478 

-0.04600 

WBm 

0.01478 

|§| 

-0.04593 

10 

1 

CF 

0.01206 

FEM 

0.01206 

mBBm 

2 

CF 

0.01019 

m 

FEM 

0.01019 

ns 

1 

CF 

0.006223 

0.4200 

-0.8655 

5 

FEM 

0.006223 

0.4200 

-0.8655 

2 

0.004975 

0.4070 

-0.09156 

H99H 

0.004972 

0.4070 

-0.09057 

Table  4.  Neutral -surface  positions  and  dimensionless  deflections 
for  cylindrical  panels  of  single-layer  (0°)  aramld-rubber 
and  polyester-rubber  under  sinusoidal  thermal  loading,  as 
determined  by  two  ‘different  methods.  (R/h  «  10.0, 
fi  »  1.0,  f  ■  0.0,  P  »  0  )* 


Aspect 

Ratio 

w 

Z 

X 

Z 

y 

a/b 

C.F. 

F.E. 

C.F. 

F.E. 

C.F. 

F.E. 

Aramld-Rubber 

0.5 

0.03016 

0.03021 

0.02250 

0.02247 

-0.07431 

-0.07308 

0.75 

0.06767 

0.06772 

0.02295 

0.02295 

-0.07458 

-0.07400 

1.0 

0.1190 

0.1190 

0.021565 

0.02158 

-0.06347 

-0.06332 

1.25 

0.1821 

0.1821 

0.01772 

0.01768 

-0.05191 

-0.05157 

1.5 

0.2545 

0.2546 

0.01124 

0.01121 

-0.04228 

-0.04204 

1.75 

0.3331 

0.3330 

0.002360 

0.02314 

-0.03473 

-0.03454 

2.0 

0.41505 

0.4149 

-0.0085575 

0.08476 

-0.02892 

-0.02876 

Polyester-Rubber 

0.5 

0.04083 

0.04068 

0.09004 

0.08958 

-0.2011 

-0.1978 

0.75 

0.08952 

0.08955 

0.08481 

0.08463 

-0.1574 

-0.1560 

1.0 

0.1527 

0.1527 

0.07445 

0.07423 

-0.1112 

-0.1109 

1.25 

0.2263 

0.2263 

0.06020 

0.06919 

-0.07815 

-0.07779 

1.5 

0.3060 

0.3059 

0.04302 

0.04300 

-0.05596 

-0.05573 

1.75 

0.3881 

0.3880 

0.02370 

0.02364 

-0.04104 

-0.04088 

2.0 

0.4695 

0.4693 

0.002851  . 

0.002742 

-0.03062 

-0.03069 

*  _ 
w  » 


Wh 


Zx  •  V 


» 


t 


zy  '  V 
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3f  ■ 


Table  5.  Neutral -surface  positions  and  dimensionless  deflections 
for  cylindrical  panel  of  two-layer  (0°/90°)  aramld-rubber 
and  polyester-rubber  under  sinusoidal  thermal  loading. 
(R/h  -  10.0,  Ti  »  1.0,  Tq  *  0.0,  PQ  »  0).* 


Aspect 

Ratio 

a/b 

m 

w 

Zx 

z 

y 

C.F. 

F.E. 

C.F. 

F.E. 

C.F. 

F.E. 

Aramid-Rubber 

1.0 

0.1212 

0.1218 

0.05189 

0.05335 

-0.05085 

-0.05133 

1.25 

0.1783 

0.1787 

0.03656 

0.03763 

-0.04511 

-0.04644 

1.5 

0.2408 

0.2410 

0.02050 

0.02084 

-0.04052 

-0.04209 

1.75 

0.3064 

0.3065 

0.003910 

0.004119 

-0.03682 

-0.03694 

2.0 

0.3726 

0.3726 

-0.01397 

-0.001406 

-0.03387 

-0.03493 

Polyester-Rubber 

1.0 

0.1829 

0;1849 

0.1486 

0.1510 

-0.09623 

-0.1006 

1.25 

0.23745 

0.2399 

0.1066 

0.1100 

-0.08727 

-0.08823 

1.5 

0.2882 

0.2905 

0.06572 

0.06857 

-0.08188 

-0.08287 

1.75 

0.33435 

0.3363 

0.02652 

0.02813 

-0.07857 

-0.08074 

2.0 

0.37525 

0.3769 

0.01112 

-0.01144 

-0.07649 

-0.07866 

w  * 
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Closed-form  and  finite-element  solutions  are  presented  for  the  thermoelastic 
behavior  of  laminated  composite  shells.  The  material  of  each  layer  Is 
assumed  to  be  thermoelastlcally  orthotropic  and  bimodular,  l.e.,  having  dif¬ 
ferent  properties  depending  upon  whether  the  fiber-direction  normal  strain  Is 
tensile  or  compressive.  The  formulations  are  based  on  the  thermoelastic 
generalization  of  Dong  and  Tso’s  laminated  shell  theory,  which  Includes 
thickness  shear  deformations.  The  finite  element  used  here  has  five  (over) 
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degrees  of  freedom  per  node  (three  displacements  and  two  bending  slopes). 
Numerical  results  are  presented  for  deflections  and  the  positions  of  the 
neutral  surfaces  associated  with  bending  along  both  coordinate  directions. 
The  closed* form  and  finite-element  results  are  found  to  be  In  good  agreement 


